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1. Introduction 

Calculation of the scattering phase shift represents an important step for expanding our un- 
derstanding of the strong interaction based on lattice QCD to dynamical aspects of hadrons. Since 
Liischer derived a finite size formula for the two-meson system on 1986 |[T]], which give us a relation 
between the phase shift and the energy eigenvalue on the finite volume, many lattice calculations 
of the scatting length and the phase shift of the two-meson systems have been carried with his for- 
mula. Recently his formula was extended to that for the Nn system by Bernard et al. by using the 
non-relativistic effective theory QCDSF collaboration calculated the phase shift of this system 
with this extended formula and study A(1232) resonance [||]. 

The extension of formula is necessary to extend our study to many systems. In the present 
work I consider a derivation of the formula for the elastic NN scattering system, where the formula 
only for spin singlet state in the non-relativistic limit, which is same as that for the two-meson 
system given by Liischer, has been known. My derivation is based only on the relativistic quantum 
field theory and any effective theories for the two-nucleon interaction are not assumed. Further the 
extension to the Nn system can be easily done as discussed latter. 

2. Wave function in infinite volume 

First we consider the wave function in the infinite volume defined by 

C^(x;k) = (0| na{x/2) pp{-x/2) |k,A„,Ap) , (2.1) 

where na(x) and p^^s) are interpolating operators of the nucleons and |k, A„, Ap) is the asymptotic 
NN state with momentum k, — k and the helicities A„, Ap . Using LSZ reduction formula, the wave 
function can be written by 

C^(x;k) =t/„/,(k,A„,Ap)e'--'^ + / I ^a/3(P,^«,^.)e''''-''^^^^^^^^g^^ , (2-2) 

where Ua^{yi,Xn,Xp) is a spinor for two free nucleons given by ?7„|3(k, A„, Ap) = MQ;(k, A„)Mj3(— k, Ap) 
with the one nucleon spinor M(k,A). r(p, i§p;k, A„, Ap) is the off-shell scattering amplitude for 
a process ?i(k,A„)/?(-k,Ap) ^ ?i(p,(§„)/7(-p,(§p). 

We can estimate ( |2.2[ ) in the region |x| > /? for the two-nucleon interaction range R, by using 
a integral formula 

/ {llyp^-i^-ie ^^P^ = ^{i-h{kx)+ni{kx))m forF(x)=0 (2.3) 

where F{x) is the inverse Fourier transformation of f{p). ji{x) is the spherical Bessel and ?i/(x) is 
the Neumann function, whose conventions agree with those in [||] as adopted in This formula 
is a extension of (A.l 1) in Appendix A in Ref. [Q] to that for arbitrary value of / and can be derived 
by similar calculations of that paper. 

From ( pj| ) we know that all values in the numerator of the integrand in ( |2.2| ) can be replaced 
by the value at on-shell p = k. The off-shell scattering amplitude r(p, (§p;k, A„, Ap) is replaced 
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by the on-shell amplitude, which can be expanded as 



T{k^p,£,„,£,p-MX^p) = 167r'xI^li.A„A,(^)-^-/'<J*(^p)<l(^'t) , (2.4) 

where Nj = {IJ + I) / {An) , A = A„ - Ap, = - <^p and ^/s = 2Vm^ + B. In (Q the helicity 
amplitude in the subspace of the total energy ^/s and the total angular momentum / is defined 



by ^^:!?,„A„A„(^) = {^nUf<^^Hk)\Xn^p). The function D\^l,iap) = (7M|exp(-/a/Jexp(-/j8/,) 
exp(— /77^)|7M') is the Wigner's D-function with the Euler angle {a,^,y) = {^p,dp,—^p) for 
momentum p = sin 0p cos 0^ , sin 0p cos 0^ , p cos dp ) . 

Using dJ) and (0), we know that the wave function ( p^ ) in the region |x| > Ris written by 



0°° (x; k) = Y^^jDmX i^k) • 07MA„ A„ ^) = ^" - 



JM 



(2.5) 
(2.6) 



where a^^}^^^^^^^{k) = ((§„(§p|/ + /f(-')/2|A„Ap) and P^^l^^^^^^^{k) = ((§„(§p|fW/2|A„Ap), which cor- 
respond to = cos 5/ •exp(/5/) and j8''' = sin 5/ •exp(/5/) for the two-meson system with the 
scattering phase shift 5/. In (2.6) the function //ma„a,,('^;'^) the wave function of two free nucle- 
ons with the total energy ^/s, the total angular momentum JM and the helicity A„Ap. Its explicit 
form is given by 



JjM^,,^i^;k)=L(V) C,,/x;^) /?(V) ^7;^,„,/x;^ 
where differential operators L{V) and ^(V) are defined by 



IR-EX 



(2.7) 



L(V) = I (g-V/O 
E + m 



E + m 



(2.8) 



with E = ^k'^ + nfi. In (2/7) the function ^ (x;/:) is 2 x 2 non-relativistic spinor defined by 
jZx„xM^k) = • {JMls\JMXnXp) , (2.9) 

Is 

Jmisi^-'k) =ji{kx)Y'/M{a,)/bi{k) , Yllt{Q.,)=Y,YUa,)^{s,pi)-C{lm;spL;JM) , (2.10) 

where the coefficient {JMls\JMXnXp) is the transformation coefficient from the helicity base to the 
orbit-spin base ((/M/5)-base) with the angular momentum / and the spin s [|6|], and C{lm;s}x;JM) 
is the Clebsch-Gordan coefficient for angular momentum state \lm) IsjJ.) and \JM). Q.^ is the 
spherical coordinate for x. bi{k) is the normalization constant of the state, which takes \/bi{k) = 
(4;r)/' • [k^ + m^) for the usual relativistic normalization {u^u = IE). <^{s,}i.) is 2 x 2 spin wave 
function for two spin 1/2 particles with total spin s^i. The function NjMi„Xp{^''>^) ( ^^ is given 
by replacing ji{kx) by ni{kx) in ( |2.7[ ). We can regard ( |2.7[ ) as a relativistic extension of the non- 
relativistic spinor ^ (x;^) to the relativistic one //ma„a,,(x;^), so that the spinor satisfies the 
Dirac equation. We use a notation /'^^Ir-ex for this relativistic extension like as (2.7) in the follow. 
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Next we rewrite ( p^ ) and ( p!q ) by the (/M/s)-base as 

r(x;k)= £QM/.(k)-&.(x;/:), (2.11) 

JMh 

^jMh (x; ^) = £ -//M/'y (x; k) ■ ajfj^^ {k) + Njmi's' (x; ^) • (k) , 



(2.12) 



with some constant Cjmis{^), where functions of the (7M/5)-base are defined by 



JjMiM'^^) = L JjMi„?i,i^;k) ■ {JMls\JMXnXp) , (2.13) 
«/v,/.(^)= L a^^l^^^^_^^^^{k)-{JMl's'\JM^„^p){JMls\JMKrXp) , (2.14) 

and NjMX„ii,{^\k) and ^f) i^{k) are similarly defined. Here we should note that JjMis{^\k) and 
A(/M/j (x; k) are not eigenstates of the orbital angular momentum and the spin with / and s. These 
functions satisfy the Dirac equation, thus the upper and the lower components have different orbital 
angular momenta. 

3. Wave function on the finite volume 

Next we consider the wave function on the finite periodic box of volume defined by 

0^p(x;^:) = (0| «„(x/2) Pui-x/l) \k) , (3.1) 

where |^) is the energy eigenstate with ^/s = iVrrP- + k^ on the finite volume. Here we assume 
the condition /? < L/2 for the two-nucleon interaction range R and the lattice size L, so that the 
boundary condition does not distort the shape of the two-nucleon interaction. In the region R < 
|x| < L, the wave function satisfies following two equations and the boundary condition. 



[ '(r- V) + fE-m] (l)^{x;k) = , ^^{x;k) \ -i{y- V) + fE 



m 



T 



0, (3.2) 



(^)^(x + nL;/t) = 0^(x;;t) ( n G Z^' ) , (3.3) 

where E = y/nP^+W. The general solution of these equations can be written by the linear combi- 
nation of the Green function defined by 

GjMU (x; k) = gZu (x; k) I , (3.4) 
GyM/.(x;^)=<(V)-^£^i-3e'P-, ^r,^(p)=y .7/^(11,), (3.5) 

where T = {p|p = (27r)/L • n , n G I?} and Q.p is the spherical coordinate for p. This Green 
function is related to that introduced in Ref . [P G/„, (x; k) by 
G^^i^{x;k) = ZmtiGim{x;k) ■ ^{s,}JL)C{lm;spL;JM). 

Using partial wave expansion of Gim(x;k) given in Ref. |jl]|, we obtain 



GjMh (x; k) = ai {k)bi (k) ■ Njmu (x; k) + ai (k) ^ Z?,/ (k) ■ JfM'i's (x; k) ■ mPj^,j, jj^j (k) , (3.6) 

J'M'l' 
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where ai{k) = (— l)'^'+'/(47r), bi{k) is the normalization constant appeared in ( 2.10| ) and 

Mi'/',7M/(^) = L M,>,n'j,„ik)-C{l'm';s^;j'M')Cilm;s^;JM) . (3.7) 

mm'll 

The function Mii^' ,im{k) in is defined by (3.34) in Ref. [[I]], which is given by 

Ml'm',lmi^) = ^l'm',l"m",lmWl»m"{(]) , q = kL/{2n) (3.8) 
l"m" 



k,n'j"rn"M = (-!)'• " [21" + 1 ) ^ • C(ZO; /"O; rO)C(/m; /"m"; /'m') , (3.9) 

WUq) = — £ ^-^^^„(n) , ^^„(n) = -YUQ-n) . (3.10) 

4. Relation between and 0^ 

In the following we restrict ourselves to the wave function for the irreducible representation of 
the rotational group on the finite volume (cubic group O), which is defined by 

#„(x;^) = (0| «(x/2) p{-x/2) \t,ra) , (4.1) 

where \k;ra) is the energy eigenstate with ^/s = l^nP- + k^ and belongs to the irreducible rep- 
resentation of O labeled by Y and a (a = 1 . . . dimP , T = {Ai ,A2,E, Ti , T2}). Projection of the 
irreducible representation of SU(2) (|/M)) to that of O (\ranJ)) is given by 

\JM) = £ \ranJ) ■ V{JM\TanJ)* , \TanJ) = ^ \JM) ■ V{JM\TanJ) , (4.2) 

Tan M 

with known coefficient V {JM;TanJ), where n is the multiplicity of the representation T. 

In the previous section the wave functions were expanded in terms of functions of the {JMls)- 



base {JjMis and Njmis)- But it is more convenient for the wave function ( |4.1| ) to expand in terms of 
functions of {ranJls)-ba.se defined by 

Jranjis{x;k) = Y,JjMh{x;k) -ViJM-XanJ) , (4.3) 

M 

with the coefficient V{JM;TanJ). 



In the region |x| > R, the wave function (4.1) can be written by the linear combination of the 



Green function and also the wave function in the infinite volume as 

#a (x; k) = Y, EvanJis (k) • GranJis (x; /:) = £ Cranjh (k) • (^YmiJh (x; k) , (4.4) 

njls njls 

with some coefficients EranJh{k) and Cranjhi^), where Gra„7/,(x;^) and ^«„y;,(x;fc) are func- 
tions of the {ranJls)-hase obtained by the transformation ( p3| ) from GjMis{x',k) defined by (| 



and ^J'j^i^{x;k) defind by ( [2.12[ ). After some calculations we obtain 



0f„(x;^) = Y,EranJhi^)[bi{k)-NranJhix;k)+ £ bi,{k) ■ Jran'j'i's{'^;k) ■ M^],,,^^ji{r;k) 

nJls \ n'J'l' , 
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= £Cra.7/.v(k)£(W/'.Kx;^)-«S,,,(r;^)+^rany/v(x;^ , (4.5) 

nJh I's' 

where the constant ai{k) are removed by redefinition of the constant Eranjisi^), and 

SrrSa'a ■Ml'],,,^j,{r;k) = ^ M^;l,,j^j{k) ■ V{J'M'-r a'n'j')V{JM-TanJ) , (4.6) 

MM' 

«/v./.v(r;^) = £a,5./, •^(■/^;ra«7)y(7M;ra«7) , (4.7) 



and pj^^J i^{r;k) is similarly defined. The diagonal property of M(r;k) in (4^6) for indices (Fa) is 

result from the invariance of M^}^^, j^^^ (k) under the rotation on the finite volume (see Ref. [[T]]). 

From (4.5), we know that coefficients of functions Jranjui^^k) and NranJisi'^',k) relate each 
other. After some calculations, we find that it is given by 

det[M(F;A:) - A(F;A:)/B(F;A:) ] = , (4.8) 

where we introduce a vector space spanned by indices {nJls) at fixed (Fa) and define linear oper- 
ators on this vector space by 

[M(r;^)] = 4s -mS,, „,,(F;/:) , [A(F;fc)]„,,,^,„,,^ = d„,,A'j ■ aSi,{r;k)/b,{k) , 

(4.9) 

and B(F;fc) is similarly defined. Equation (1^) is a finite size formula for the elastic NN scattering 
system, which gives us a relation between the energy eigenvalue on the finite volume and the 
quantity of the elastic scattering A/B. 

5. Finite size formula for A^A^ system 

In this section we show the explicit matrix form of the finite size formula for the A^A'^ system 



(|4.8|). S'-matrix at fixed J forms a 4 x 4 matrix. This matrix is reduced to sub-matrices by the 
eigenvalue of the global symmetry : the parity P and the particle exchange R (= (—1)^ with the 
iso-spin /) as 

S^^^ = (2x2matnx;P = {-iy-'^,R = {-\y-'^). (5.1) 
+ M X 1 matrix ; P = {-ly , R = {-\y 
+ ( 1 X 1 matrix ; 7^ = (-1)^ ,R = {-iy-^ 



A(F;^) and B(F;^) in the finite size formula (O) also take same form 



We note that the basis of the partial wave expansion Jranjisi'^^k) and Nrmjisi^'^k) in ( p3| ) are 
eigenstates of the parity and the particle exchange with /" = (—!)' and R= (— 1)' • (— 1)*^'. Thus 
the wave function for the state with R = —P, only functions with = appear in the partial wave 
expansion. For the state with R = P, only functions with = 1 appear. The mixing between s = 
and 5 = 1 is forbidden by the symmetry of the parity and the particle exchange (iso-spin). Therefore 
we can separately obtain the finite size formula for P = —R {s = 0) and P = R (s = 1). 

In the case of R= —P (s = 0), the components of the matrix M(F;/c) in the finite size formula 
( P^ ) are given by 

^nii'M (r; k) = 5jn' 5j, ■ £ MfM',jM (k) ■ V (Fa«/; /M') V (TanJ; JM) . (5.2) 

MM' 
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This is the same matrix as that appeared in the finite size formula for the two-meson system. 
Further, af^i^{k) = djiSjv ■ ai{k) and ^^^^^^{k) = SjiSjv ■ pi{k) with the diagonal components ai{k) 
and pi{k) also take same matrix form as that for the two-meson system. Thus the finite size formula 
for the A^A'^ system with P = —R {s = 0) is same as that for the two-meson system in Ref. 

In the case of R = P (s = 1), the matrix M{r;k) and A{r;k) /B(r;/c) have complicated struc- 
ture. In the following we show explicit matrix form of the finite size formula for some channels as 
example. We neglect contributions of / > 5. In this case the multiplicity n is 1 for all irreducible 



representations F, thus we omit the index n in the formula ( |4.8| ) for simplicity, as 



det[M(F;A:) - A(F;A:)/B(F;A:) ] = 



(5.3) 



[M(r;^)],,,,,,=M,i;>,,„,;(r;^) , [A{r;k)]j,,,^j^ = 5rj-alfJ,^{r;k)/bi,{k) , (5.4) 

where n = n' = I, s = s' = I and the matrix B(F;/c) is similarly defined. 

The first example is the deuteron state. We have to consider the NN state with the total angular 
momentum / = 1 and the parity P = +1, which corresponds to ^Si and ^Di states in the non- 
relativistic limit. The J = I state belongs to the irreducible representation of the cubic group 
r = Ti, thus we consider the finite size formula for F = Ti for the study of the deuteron. The other 
angular momentum states also belong to Ti as Ti = 1 + 3 + 4 up to 7 > 5 and the finite size formula 
includes contributions from all these states. For each values of J, possible values of I are given by 



1 = 0,2 for / = 1 
1 = 2,4 for / = 3 
1=4 f or / = 4 



(5.5) 



from the parity conservation and the theory of addition of the angular momentum. Thus matrices 



M(F;/t) and A(F;it) /B(F;/t) in the finite size formula (5.3) take : 



M 



/ MiOJO Mio,12 MiQ32 Mio,34 MiOM \ 
MuAQ M\2A2 M[232 ^12,34 ^12,44 
■^32,10 -^32,12 -^32,32 ■'^32,34 -^32,44 
Af34,10 ^34,12 ^34.32 ^34,34 -^34.44 

VM44^10 M44J2 M44.32 -^44,34 ^44.44 / 



A/B: 



















J=3 














J=4 J 



, (5.6) 



where the boxes in the matrix A /B which enclose the values of / refer to 2 x 2 or 1 x 1 matrices 
expanded by the possible values of /. In (5^) components of the matrix M are denoted by Mjqiji = 
\^(jr;k)]jii'ji and are given by 



Mio,io = Woo 
Mi2,io = 
M32,10 = 
M34,10 = -2W40 

M44,10 = f VVW. 



40 



Mi2,12 
^32,12 
^34,12 
M44J2 



Woo 

-f\/6W4o 

^W40 



M32,32 
^34,32 
^44,32 



Woo + fW40 



§^3^40 + 1^3^60 



18 
77 



21W40 + |^V2W60 



11 



M34,34 =Woo + f^W4o + ffW60 

M44,34 = - If VlWeo M^.u = Woo + ^ W40 + ^Weo + -^Wgo , (5.7) 
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where Mjip ji = Mji^n' and the function W/m is defined by (3.10). 

Finally we consider the state with same J but opposite parity to the deuteron, ie. J = I, P = —I. 
We also consider the representation F = Ti . Possible values of / for each J are different from those 
of the deuteron case as 

for 7 = 1 

for 7 = 3 . (5.8) 
for 7 = 4 



1 = 1 
1 = 3 
1 = 3,5 

Thus matrices in the finite size formula take different forms 

/ Mn,n Mn,33 MnA3 Mn,45\ 

M33J1 M33^33 M33^43 M33^45 
M43,ll M43,33 M43,43 M43,45 
\M45a1 M45,33 M45,43 M45,45 / 

^33,11 

^43,11 = - 7 V^iUVl'40 
■M45,ll = 7 V'+ZV1'40 



M: 



A/B 



J=l 













J=3 













J=4 



Woo 

^jVl4W4o 
-iV2T0W4i 



Woo + ^iW4o - fiWeo 



^33,33 

M43.33 = -^^^40 - i V15W60 
M45,33 = fl^W40+'iV3W60 



43,43 — W'oo + JjW^o ■ 33 , 

M45,43 = V5W40 - T^VSWeo - ^VSW^o 



45,45 



Woo + ^W40 - ^^60 - 7x1^80 



(5.9) 



(5.10) 



6. Summary 

The finite size formula for the elastic A^A^ scattering system is derived from the relativistic 
quantum field theory. The extension to other two-baryon system as the A^A system is trivial. Finally 
I give a comment for the Nn system. The formulation of this paper for the A^A'^ system is also valid 
for the system with the general value of the spin s. Thus the finite size formula for the Nn system 
can be easily obtained from that for the A^A'^ system (18) by set = 1 /2. In calculations of matrices 



M(F; k), A(F; k) and B(F; k) in (g^, we change the coefficient y(7M; Fa?i7) in (gj2|) by that for the 
double covered cubic group (^O) to deal with the half integer value of the total angular momentum 
J as discussed in Ref. [^. I confirmed that my results are consistent with those obtained from the 
non-relativistic effective theory by Bernard et al. [^. 

This work is supported in part by Grants-in-Aid of the Ministry of Education (No.20540248). 
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